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RESUMEN 

En este trabajo estudiamos las propiedades de la emision en radio de un conjunto de galaxias HII, con el 
proposito de detectar regiones de formacion estelar jovenes, en las cuales las primeras explosiones de supernova 
estan teniendo lugar. Las distribuciones espectrales de energfa observadas tienen tres comportamientos difer- 
entes: 1) hay galaxias con una distribucion de energfa del tipo sincrotron, 2) galaxias con espectro termico y 
3) galaxias que presentan absorcion libre-libre a longitudes de onda larga. Este ultimo caso es un indicador 
de la presencia de un cumulo estelar masivo que se cncuentra aun dentro de la nube molecular progcnitora. 
Comparando las tasas de formacion estelar (SFR) determinadas a partir de lineas de recombinacion con las 
tasas determinadas a partir de las observaciones de radio encontramos que el SFR(Hq;) es en promedio 5 voces 
mayor que el SFR(1.4 GHz). Estos resultados sugieren que la emision dc estas galaxias esta dominada por un 
episodio de formacion estelar reciente y masivo en el cual las primeras supernovas estan empezando a explotar. 
Concluimos que el deficit de emision sincrotron en aquellas galaxias con los mayores anchos equivalentes de 
H/3 solamente puede ser explicado si estos brotes de formacion estelar son extremadamente jovenes con edades 
menores que 3.5 millones de afios de forma que aun no lian cxplotado las primeras supernovas de tipo H. 

ABSTRACT 

We studied the radio properties of very young massive regions of star formation in HH galaxies, with the aim 
of detecting episodes of recent star formation in an early phase of evolution where the first supernovae start to 
appear. The observed radio spectral energy distribution (SED) covers a behaviour range; 1) there are galaxies 
where the SED is characterized by a synchrotron-type slope, 2) galaxies with a thermal slope, and 3) galaxies 
with possible free-free absorption at long wavelengths. The latter SED represents a signature of massive star 
clusters that are still well inside the progenitor molecular cloud. Based on the comparison of the star formation 
rates (SFR) determined from the recombination lines and those determined from the radio emission we find 
that SFR(Ha) is on average five times higher than SFR(1.4 GHz). These results suggest that the emission of 
these galaxies is dominated by a recent and massive star formation event in which the first supernovae (SN) 
just started to explode. We conclude that the systematic lack of synchrotron emission in those systems with 
the largest equivalent width of H/3 can only be explained if those are young starbursts of less than 3.5Myr of 
age, i.e. before the first type H SNe emerge. 
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1. INTRODUCTION 

Thanks to their proximity and the posibillity of 
observing them in almost all the electromagnetic 
range, nearby galaxies are unique laboratories to 
study galaxy evolution. This contribution focus on 
the study of early stage of stelar evolution when most 
of the new formed stars remain deep inside the pro- 
genitor molecular cloud. In fact, during the early 
stages of evolution, very young star clusters should 
either lack or have a deficit of synchrotron emission, 
being dominated by free-free radio emission. Given 
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that the synchrotron emission observed in starbursts 
is due to supernova activity, it is a direct tracer of 
the end of the evolution of massive stars (M > 8 Mq). 
In the case of a coeval star formation burst, the first 
type II SNe appear after only 3.5 Myr and last imtil 
the 8M0 stars explode at around 40 Myr. Conse- 
quently, the lack of synchrotron emission in presently 
starforming regions is a good indicator of very young 
clusters. HII galaxies are low mass objects whose 
emission and thus most observables are dominated 
by a massive burst of star formation. Their optical 
spectrum is identical to that observed in giant HII 
regions like 30-Doradus in the LMC. These proper- 
ties make these galaxies ideal targets for a study of 
integrated properties in search for features related 
to massive star clusters. Selected sources can then 
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be studied in more detail with higher resolution ob- 
servations. In this sense HII galaxies are considered 
"young" . In fact they are probably the youngest stel- 
lar systems that can be studied in any detail. The 
youth scenario is supported not only by the strong 
emission lines from ionized gas, but also by their un- 
derlying stellar continuum properties, i.e. no stellar 
absorptions are detected in those HII galaxies with 
the strongest emission hues (EW(H/3) >150A) and 
only weak hydrogen and helium and some very weak 
metal lines are detected on those more evolved HII 
galaxies (EW(H/3) < 50A). The lack of detection of 
stellar absorptions in the extreme HII galaxies is due 
to the fact that the optical spectrum of very young 
clusters is dominated by the light of massive blue su- 
pcrgiants close to the main sequence turn-off. The 
spectrum of these stars shows narrow and relatively 
weak absorptions of hydrogen and helium plus ex- 
tremely weak metal lines. In the case of HII galaxies 
(and also HII regions) the weak and narrow stellar 
absorptions are filled with strong emission from the 
ionized gas of the associated HII region making it 
impossible to detect directly their presence. 

Here we present the results of a study of inte- 
grated radio spectral energy distributions (SEDs) of 
nascent starburst candidates. We find several galax- 
ies with flat spectra, typical of thermal emission, or 
even with inverted spectra, indicating the presence 
of heavily embedded star clusters. With the data at 
hand we compare the radio and emission line proper- 
ties of these galaxies and propose a qualitative evo- 
lutionary model. 

The long term objetives of this project includes 
not only to find an d study the youngest bursts i n the 
nearby Universe (jRosa-Gonzalez et al.l l2007al ) but 
also to understand the transition between the dusty 
obscured {z between 2 and 3) to the optically bright 
~1) Universe and the potential use of H II galaxies 
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as standard candles ( Melnick et al. hood) 



2. FIRST EVIDENCE OF DEALING WITH 
NASCENT GALAXIES 

Our sample selection was based on the expec- 
tation that in young starforming galaxies their most 
massive stars did not have enough time to evolve and 
explode as supernovae. In a situation like this, these 
galaxies should have very little synchrotron emission 
([Bressan et al.l 120021 ). Consequently, their radio lu- 
minosities should be smaller than what one would 
expect to find based on the SFR measured from other 
indicators, like Ha. 

For the current study we sele cted a sample 
from the catalog of HII galaxies by iTerlevich et al 
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Fig. L Evolution of the EW(H/3) for an instantaneous 
burst, solid line, and for continuous star formation, 
dashed line. 



(|l99lf ). Based on the reddening corrected Ha fluxes 
of the galaxies in this catalog, we estimate their 
St ar Formation Rates (SFR) u sing the relation given 
bv lRosa-Gonzalez et all (|2002l) . 



SFR[Moyear"^] = 1.1 x 10""*' x L(Ha) (1) 

where the L(Ha) is given in erg s~^. These SFRs 
were then converted into radio 1.4 GHz fl uxes, using 
the relation given bv lSchmitt et al. ( 2006[ ). Compar- 
ing the predicted 1.4 GHz fluxes with values obtained 
from the NRAO/VLA Sky Survey (NVSS) catalog, 
we selected those galaxies which were not detected 
at the NVSS 5ct limit (2.5 mJy), or have an observed 
flux smaller than the predicted one. This deflcit of 
1.4GHz emission can be quantiflcd in terms of the 
ratio between the observed 1.4 GHz flux and the ex- 
pected flux based on the SFR calculated using the 
optical emission lines. For simplicity we denoted this 
ratio as the d-parameter. 

The EW(H/3) has long been us ed to es t imate 
the age of a stellar bur st (e.g. iDottoril Il98ll 
Luridiana fc Peimbertll200l[) . However, the EW(H/3) 
indicator which consists of the relation between the 
continuum flux which depends on the whole star for- 
mation history of the galaxy and the strength of the 
emission line which depends on the recent (^ 5 x 10^ 
years) star formation activity can be signiflcantly 
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lower than the age of the current burst. This prob- 
lem has been noticed i n the analysis of the most 
extreme HII galaxies bv iTerlevich et alJ (|2004l ). To 
illustrate this problem, we plot in Figure[T]the evolu- 
tion of the EW(H/3) for the case of an instantaneous 
burst and for continuous star formation. Both mod- 
els were calculated using the SB99 code for the case 
of a Salpeter initial mass function and masses be- 
tween 0.1 and IOOMq. In both modes the EW(H/3) 
decreases with time but due to the formation of new 
massive stars in the continuous mode, the EW(H/3) 
remains larger as time increases. These star form- 
ing histories mark the two extrema and the age of a 
galaxy with a given EW(H/3) must lie between the 
limits shown in Figure [TJ In the case of the exis- 
tence of an old population, the observed EW(H/3) is 
reduced even more - in comparison with the contin- 
uous case - due to the integrated light of the galaxy 
and the absence of massive stars responsible for the 
recombination lines. Interesting enough is the fact 
that there are no galaxies known to have EW(H/3) 
greater than about 350A. 

Figure [2] shows the EW(H/3) against the d- 
parameter expressed in percentage. Notice that most 
of our galaxies lie in a region where the observed 
flux is less than 50% of the expected value, showing 
that indeed our sample has a deficit of synchrotron 
radiation. Five galaxies (Tol 116-325, MRK1315, 
Tol 1303-281, Tol 1304-386 and Tol 1358-328) arc ex- 
treme cases with the highest EW(H/3) and just upper 
limits in the observed 1.4 GHz fluxes. Notice that, 
based on the EW(H/3), these galaxies have ages of 
less than 10^ years even in the case of continuous 
star formation (see Figure [T]). 

In order to investigate how our galaxies compare 
to normal, quiescent star forming galaxies in the 
d— parameter us.EW(H/3) diagram, we obtained ra- 
dio 1.4 GHz fl uxes (from NVSS) fo r a sample of suc h 
galaxies from Jansen et al. I (|2000^ : lHo etall (Il997l) . 
We can see in Figure [H that the normal galaxies 
have small EW(H/3) and large d values, occupying 
a different region relative to the galaxies in our sam- 
ple, again consistent with the interpretation that our 
galaxies are young and have a deficit of synchrotron 
radiation. A particularly interesting result from this 
Figure is the fact that a significant portion of the nor- 
mal galaxies, those with the smaller EW(H/3), have 
d >100%. Although this result seems contradictory, 
since it implies a higher radio SFR than what is de- 
rived from Ha, it can be explained in view of the 
lifetime of these two indicators. For example, in the 
case of a single burst, the Ha emission lasts only 
10 Myr, imtil all the ionizing photons die. The ra- 
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Fig. 2. H/3 equivalent width against the d parameter de- 
fined as the ratio between the observed and the expected 
1.4 GHz flux (see text for further details). Circles are 
the sample HII galaxies and open symbols represent the 
control sample from Jansen et al. (2000) (triangles) and 
from Ho et al. (1997) (squares). 

dio emission on the other hand will last much longer 
than that, because of the longer synchrotron emis- 
sion life time and the fact that the stars that explode 
as SN last for ~40 Myr. 

3. RADIO SPECTRAL ENERGY 
DISTRIBUTION 

The galaxies were observed with the VLA at 
4.9 GHz (6 cm) and 8.4 GHz (3.5 cm) and we comple- 
mented the data set by including the 1.4 GHz data 
from the NVSS catalog. In Figure [3] we present the 
radio SEDs of three selected galaxies. Analyzing the 
spectral slopes of the galaxies detected at 1.4 GHz 
we find that they have spectral indices between those 
expected for thermal emission (a = —0.1), and those 
characteristic of a source dominated by syn chrotron 
radiation (typical value is around a = — 0.8. [c"ondon 
1992. 

In Figure U] we present the distribution of our 
galaxies in the F(1.4 GHz) vs. F(8.4 GHz) dia- 
gram. This Figure also presents two lines indicating 
the location of thermal (a = —0.1) and synchrotron 
(a = —0.8) emission. In the case of more evolved star 
forming galaxies, where we find a mixture of these 
two components, the 1.4 GHz emission is dominated 
by synchrotron emission, while at 8.4 GHz both syn- 
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Fig. 3. Radio SED of three representative galaxies. The SED of UM533 shows the effect of the free-free absorption 
which have reduced significatively the flux at 1.4 GHz. Mrkl318 shows a typical, almost flat, spectrum of a thermal 
source. UM488 shows a steep spectrum consistent with being dominated by synchrotron emission. 



chrotron and thermal emission produce similar con- 
tributions. Most of the galaxies in our sample lie in 
the region between these two lines, indicating that 
they are young objects still dominated by thermal 
emission. 

It is obvious that Figure 5] represents an oversim- 
plification, as the history of star formation is much 
more complicated and cannot be represented by ei- 
ther continuo us or a single burst or even by a couple 
of bursts (e.g. lXerlevich et al.ll2004[ l. Keeping this in 
mind we overplot two qualitative evolutionary tracks 
to the points in Figure^ In the first case (solid thick 
line) we assumed that the galaxies start with only 
thermal emission from the OB associations, which 
are heavily absorbed at both 1.4 GHz and 8.4 GHz 
and can only be detected at the highest frequencies. 
In this phase the dense giant HII regions arc char- 
acterized by optically thick free-free emission, which 
is observed as an inverted spectrum (thermal) radio 
source and it is commonly named ultradcnsc HII re- 
gion. Galaxies for which we only have upper limits 
at 1.4 GHz, seem to be consistent with that descrip- 
tion. 

Due to the effect of stellar winds, with typical 
velocities of 10'^ km s~^ and mass losses of about 
3 x10~^Mq year~^, the absorbing material starts 
to blow away, the optical depth is reduced and a 
galaxy moves in this diagram towards the direction 
of UM533 (4.2;2.6, in Figured]). Once all the ab- 
sorbing material is blown away, a galaxy moves up 
in this diagram, until it reaches the line where the ra- 
dio emission is dominated by thermal radiation. As 
the galaxy continues to evolve, SNe explode produc- 
ing copious amounts of synchrotron radiation and 
moving the galaxies to the top right portion of this 
diagram wh ere the relation b etween radio and FIR 
is obtained (jYun et al.l 120011 ) . In this scenario we 



expect young galaxies to have a > —0.1, with the 
younger objects having the larger slopes due to free- 
free absorption, evolving towards a = —0.8. 

In the second scenario we start with some syn- 
chrotron emission due to a previous generation of 
supernovae. If a new burst occurs a new generation 
of stars is produced, increasing the amount of ion- 
izing photons, and moving the galaxy to the right 
towards the thermal relation (see dashed thick line 
in Figurc[4]). When the next generation of supernova 
explodes the galaxy starts to move up and right to 
the position of the non-thermal relation (dot-dashed 
line). 

4. COMBINING OPTICAL AND RADIO 
MEASUREMENTS 

In this section we combine the radio and optical 
data described in previous sections in order to test 
the evolutionary models outlined in section [3] 

After an episode of star formation, the first SN 
explosion and consequently the first contribution to 
the synchrotron emission appears at about 3.5x10^ 
years, which corresponds to an EW(H/?) of about 
120 A for an instantaneous burst or 310 A if star for- 
mation is continuous (Figure [T]). Therefore objects 
with high EW(H/3) and flat slopes are candidates 
to be galaxies dominated by very recent starburst 
episodes. 

Massive stars within a young burst ionize the 
surrounding media producing free electrons and the 
presence of recombination lines. The intensity of the 
H/3 emission line is related to the free-free radiat ion 
observed at radio wavelengths by (|Condonlll992l ). 



H/3 (ergs-icni-2) = 2.8 x 10-"(Te/104)-°-52j,° iF, 

(2) 
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F(0.4 GHz), mJy 

Fig. 4. Distribution of F(1.4 GHz) as a function of 
F(8.4 GHz) for the selected galaxies. The dashed thin 
line shows the expected position of galaxies dominated by 
thermal emission, while the dot-dashed line shows where 
galaxies dominated by non-thermal emission are located. 
Superimposed on this plot are solid-thick and dashed- 
thick lines indicating qualitative evolutionary tracks de- 
scribed in the text. 

where Te is the electron temperature of the plasma, 
V is the frequency of the observation (expressed in 
GHz) and ¥^ is the flux in mJy at the given fre- 
quency. For galaxies dominated by synchrotron ra- 
diation, the thermal emission is about 10 percent 
of the synchrotron radiation at 1.4GHz. However, 
due to the different behaviour and slopes of thermal 
and non-thermal emission, the thermal emission at 
8.4GHz (3.5cm) could be about 50 percent of the 
synchrotron one or even more for those galaxies that 
are synchrotron deficient. Figure[5] (top panel) shows 
the extinction corrected H/3 flux against the flux at 
8.4GHz. Most of the galaxies have fluxes which are 
consistent with thermal emission from plasmas with 
temperatures between lO^K (solid line) and 5xlO'*K 
(dashed line) . However some galaxies have an excess 
of radio emission due to the presence of synchrotron 
emission. 

At 20 cm (1.4GHz: bottom panel in Figure[5]) the 
flux is clearly contaminated by the presence of syn- 
chrotron emission. As in the top panel the solid and 
dashed lines represent the thermal emission based on 
Equation [2l All the galaxies detected at 20 cm have 



fluxes above the thermal limits, however the galaxies 
not detected at 20 cm (the upper limits) are close to 
the thermal region. 

Due to the presence of massive stars in normal 
and starburst galaxies, there is a strong relation be- 
tween the observed FIR and the flux at 1.4 GHz. 
This correlation, which covers several orders of mag- 
nitude in luminosities, is explained by the presence 
of massive stars which heat the dust producing the 
observed FIR radiation and, as a consequence of the 
short life time of these stars, they explode as SN 
produ cing the synchro tron emission observed at 1.4 



GHz (|Yun et al.l 120011 ). Combining this strong cor- 



relation with the relation between the FIR luminosi- 
ties and the current SFR from Kennicutt (1998), Yun 
and collaborators proposed a robust relation between 
the luminosities at 1.4 GHz and the SFR, 



SFR[MQyear-^] 5.9 x 10"^-* x Li.4ghz (3) 

where the luminosities are expressed in units of erg 
s^^ Hz^^. Notice that all SFR estimates have un- 
certainties, because each estimator traces different 
populations and depends on several assumptions, 
e.g. the initial mass function, or on the light at- 
tenuation which is a strong function of wavelength. 

Comparing the SFR estimated by the 1.4 GHz 
luminosities with those given by the recombination 
lines we find that the optical SFR is on average 4 
times higher than the SFR given by the radio lumi- 
nosities. In the bottom panel of Figure [5] we plot a 
dot-dashed line where the optical SFR is equal to the 
radio SFR based on the Yun et al. (2001) relation. 
Only three galaxies are close to the dot-dashed line 
showing that most of the observed starbursts still are 
in the early stages of evolution, when not too many 
stars have become SN, and as a result the observed 
radio SED is dominated by thermal emission. The 
most extreme cases are those galaxies with non de- 
tection at 20 cm. The cases with clear signature of 
free-free absorption and those dominated exclusively 
by thermal emission are similar to regions detected 
in NGC625, and in He2-10. 

5. DISCUSSION AND CONCLUSIONS 

We combined optical spectroscopy with radio 
continuum observations to detect systematically ex- 
tremely young systems, before the first Type II SNe 
explode. The lack of these energetic events produce 
a deficit in the observed radio fluxes that was quan- 
tified introducing the d-parameter. 

To fully explain the observed trends, a highly 
synchronized star formation is required, i.e. the mas- 
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Fig. 5. Relation between the extinction corrected H/3 
flux, the flux at 8.4 GHz (top panel) and the flux at 1.4 
GHz (bottom panel). In both panels the solid line is 
the relation between the strength of the H/3 line and the 
corresponding radio flux for the case of thermal emission 
coming from a plasma with a temperature of lO'^K. The 
dashed lines correspond to a temperature of SxlO^K. In 
the botto m panel we draw (dot-dashed-line) the relation 
given by lYun et af] |20o3) which includes synchrotron 
emission. In both panels, the solid symbols are galaxies 
with a > —0.5. Non detections are represented by the 
arrows (upper limits). 



must be included in the estimation of photometric 
redshifts of distant obscured galaxies that will be 
discover with the near future facilities (e.g. ALMA 
and LMT/GTM). In fact, high redshift (sub)mm- 
galaxies with extremely weak optical counterparts 
could be the analogous of the galaxies discussed in 
this talk (with a stellar mass ~100 times higher !!!). 

Because the high energy (2 - 10 keV) emis- 
sion of star-forming galaxies is mainly due to 
HMXBs, for a given SFR, the lack of X-ray emis- 
sion could be a sign of a young system in which 
the pop ulation of compact objects is not fully de- 
veloped ( Rosa- Gonzalez et al. 2007bl) . 

An extensive version of this work that includes 
details of the VLA observations, together with fur- 
ther discussions which incorporate among others the 
analysis of the FI R emission and the q-paramete r has 
been presented in iRosa-Gonzalez et al.l (|2007al) . 
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